Memory effect on the formation of drying cracks 



in 
o 
o 

(N 
> 

o 

00 



C/2 



Michio OtsukO 

Department of Pure and Applied Sciences, University of Tokyo, 

(Dated: February 2, 2008) 



Komaba, Tokyo 153-8902, Japan 



We propose a model for the formation of drying cracks in a viscoplastic material. In this model, we 
observe that when an external force is applied to a viscoplastic material before drying, the material 
memorizes the effect of the force as a plastic deformation. The formation of the drying cracks is 
influenced by this plastic deformation. This outcome clarifies the result of a recent experiments 
which demonstrated that a drying fracture pattern on a powder-water mixture depends on the 
manner in which an external force is applied before drying. We analytically express the position of 
the first crack as a function of the strength of an external force applied before drying. From the 
expression, we predict that there exists a threshold on the strength of the force. When the force 
applied is smaller than the threshold, the first crack is formed at the center of the mixture; however, 
when the force applied exceeds the threshold, the position of the first crack deviates from the center. 
The extent of the deviation increases as a linear function of the difference between the strength of 
the force and the threshold. 

PACS numbers: 46.50.+a, 83.60.La, 46.35.+Z 
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I. INTRODUCTION 

Cracks are observed on various materials such as rocks, 
tectonic plates, and paintings. These cracks are fasci- 
nating and has been studied by many researchers. The 
study of drying cracks in a powder-water mixture was 
also included in these researchers. In one example, it 
was noticed that when a layer of a powder- water mixture 
is dried in a container, it shrinks and cracks are formed 
on it Q . The cracks extend from the surface of the mix- 
ture to the bottom and propagate horizontally along a 
line. As a result, a two-dimensional fracture pattern is 
formed on the surface of the mixture. 

When we gently pour a powder-water mixture into 
a container and leave it undisturbed during the drying 
process, a random, isotropic fracture pattern is formed. 
However, Nakahara and Matsuo reported that when an 
external force is applied to the mixture before drying, 
the fracture pattern changes depending on the manner 
in which the force is applied For example, when 
the mixture is vibrated in one direction before drying, 
cracks that are perpendicular to the direction of the vi- 
bration emerge first. Finally, a lamellar fracture pattern 
is formed. It takes more than 3 days for the cracks to be 
formed after the vibration. It is quite surprising that the 
effect of applying the force remains for such a long time. 

The experimental result attained by Nakahara and 
Matsuo represents that a fracture pattern is controlled 
by the memory of an external force applied before drying. 
Similar memory effects that a response can be controlled 
by the memory of an operation have been observed in 
other materials, such as sand piles J3, micro-gel pastes 
0, and rubbers j^. By recalling that these memory ef- 
fects have been studied from a rheological point of view. 
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we conjecture that the rheological property of a mixture 
plays an important role in the memory effect on the for- 
mation of drying cracks. 

Among the rheological properties of the mixture, the 
most conspicuous one may be plasticity. Hence, we study 
the role of plasticity in the memory effect on the forma- 
tion of drying cracks. First, in Sec. ^1 we propose a 
model of the formation of drying cracks in a viscoplastic 
material. In Sec. IIIIL we find that a viscoplastic ma- 
terial memorizes the effect of an external force before 
drying as a plastic deformation by calculating the model 
numerically. By the influence of this plastic deforma- 
tion, when an external force is applied to the material, 
a crack perpendicular to the force emerges earlier than 
when no force is applied. Based on this result, we conjec- 
ture that the perpendicular cracks emerge first by a plas- 
tic deformation. Furthermore, we express the position of 
the first perpendicular crack in terms of measurable ma- 
terial properties in Sec. IIVI This result can be used to 
test our conjecture that the memory effect is caused by 
a plastic deformation. Section [3 is devoted to the sum- 
mary and discussion. Technical details are summarized 
in Appendix 1X1 



II. MODEL 

We propose the model of a viscoplastic material in a 
similar way as that demonstrated by Ooshida and Seki- 
moto ^ . We consider a viscoplastic material of thickness 
H and width 2L in a container, as shown in Fig. ^ The 
coordinate system (x, z) is assumed such that the center 
of the container is at a; = and the bottom is at z = 0. 
For the mathematical simplicity, we restrict our attention 
to plane strain deformations of the viscoplastic material 
and we consider only a displacement u{x,z,t) in the x 
direction and a plastic strain s(x, z, t), which express the 
occurrence of a plastic deformation. We assume that the 
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FIG. 1: The illustration of viscoplastic material in a container. 

time evolutions of u{x, z, t) and s{x, z, t) are described as 
du{x,z,t) dcrxxix,z,t) daxz{x,z,t) 



7- 



dt 



dx 



+ 



dz 



a{tl (1) 



B 



ds{x, z, t) 
Ft 







i\<J..\-<JY{t)) 



Wxz\ < <JY{t) 
<^xz (2) 



otherwise, 



where 7 is a coefficient of a viscosity, a{t) is an external 
force, cry (i) is a yield stress, axx{x^ 2, t) is a normal stress, 
axz{x, z,t) is a shear stress and i? is a coefficient which 
determines the speed of the plastic strain. <Txx{x, z, t) and 
axz{x,z,t) are determined by constitutive equations: 

axx{x,z,t) = {\ + 2^l)(^^^^^^^ + c{t)^, (3) 

, . / du(x, z,t) \ 
axz{x,z,t) = — s[x,z,t)\, (4) 

where A and /i are lame coefficients, and c{t) is a refer- 
ence strain rate. An increase in the reference strain rate 
represents an effect of shrinking by drying 0, ■ 

We consider the case where the functional forms of 
a{t), ayit), and c{t) are given as 



a(t) 



_ J sin(7rt/ri) t < Ti 
t>Ti, 



00 t > T2, 



c{t) 



t<T3 
h{t-T^) t>T3, 



(5) 



(6) 



(7) 



where om, o'yoj a-nd b represent the maximum value of 
an external force, the yield stress before drying, and the 
speed of drying, respectively. 

In order to specify the model completely, we assume 
the boundary conditions. At the region where the mate- 
rial comes into contact with the container, the displace- 
ment is set to zero. In contrast, at the region where 
there is no contact between the material and the con- 
tainer, the stresses applied to the free surface become 



zero. Here, one peculiar phenomenon arises: when the 
material is dried, it peels off from the walls of the con- 
tainer. This phenomenon implies that the boundary con- 
ditions change when the peeling occurs. Hence, we as- 
sume boundary conditions 



u{±L,z,t) = 0, 
u{x,0,t) = 0, 
(^xz{x,H,t) = 0, 



(8) 



for < i < T2, and 



axxi±L,z,t) = 0, 
u{x,0,t) = 0, 
(Txz{x,H,t) = 0, 



(9) 



for t > T2, where we consider that the peeling occurs at 
t = T2. In Tabled we summarize the functional forms of 
the parameters and the boundary conditions at the walls 
of the container when Ti < T2 < T3. 

Finally, we assume the condition of a crack forma- 
tion. To the best of our knowledge, the condition has 
not yet been completely understood; however, two con- 
ditions have been used in previous works. These are the 
critical stress condition and the Griffith criterion. Un- 
der the critical stress condition, a crack is formed when 
the stress exceeds a material constant. This has been 
used in several numerical models because of the tech- 
nical advantage of the local condition. In contrast, under 
the Griffith criterion, a crack is formed when the energy 
released during the formation of a crack exceeds the in- 
crease of the surface energy. This condition has been 
used by Komatsu and Sasa in their theory Q . Although 
we cannot determine which condition is more efficient, we 
employ the critical stress condition for the simplicity of 
the treatment. Concretely, we define an average normal 
stress at the position x as 

{axx{x,t)) ^ ^ J dz axx{x,z,t). (10) 

Then, the condition of a crack formation is given as 

{a-xxixc,tc)) ^ (Jb~^ (Txx{xc,z,t) = for i > i^, (H) 

where Xc is the position of a crack, and tc is the time 
when the crack is formed. We only consider cracks that 
are perpendicular to the x axis. 

To summarize, our model consists of Eqs. and Q 
with parameters given by Eqs. ((SJ, JHl, and |(7J under 
the boundary conditions JHl and 0, and the condition 
of a crack formation Hll|l . 



III. A QUALITATIVE COMPREHENSION OF 
THE ORIGIN OF THE MEMORY EFFECT 

In order to understand how a viscoplastic material 
memorizes the influence of an external force that is ap- 
plied before drying, we numerically calculate Eqs. 
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TABLE I: Boundary conditions (B.C.) and parameters 





t 


a 


c 


ay 


B.C. at walls 


Applying an external force 


< t < Ti 


aM sm{nt /Ti) 





CTYO 


u{±L,z,t) = 


Relaxation 


Ti <t<T2 








cyo 


u{±L,z,t) = 




T2<t<T3 








oo 


(Txx{±L,Z,t) = 


Drying 


n < t 





b(t - n) 


oo 


axx{±L, z,t) =0 



and ||2J) with the control parameter q;j\/. In this nu- 
merical calculation, the initial conditions are given as 
u{x,z,0) = s{x,z,0) = 0, and the parameter values are 
set as = 1.0, L = 10.0, A = 1.0, = 0.1, aya = 0.05, 
at = 0.01, 7 = 1.0, B = 1.0, Ti = 30, T2 = 60, T3 = 90, 
and b = 0.00019. In the following we report all quantities 
in reduced units, i.e., length in units of H, stress in units 
of A , time in units of (jH'^/X) and other quantities in 
units of the combinations of these units. 

When no external force is applied to the material be- 
fore drying, (that is, aM = 0) stresses become zero at 
t = T3 (just before starting drying process). Even when 
an external force is applied to a material, provided aM is 
not so large that the shear stress exceeds the yield stress, 
the stresses remain zero at t = T3. However, when aM 
is sufficiently large, the stresses have a non-zero value at 
t ^ T3, by a plastic deformation. 

As an example, we show the results in the case aM = 
0.08. When an external force a(t) increases in time from 
0, the material displaces to the direction of shown in 
Fig. [3 Due to the influence of the boundary conditions 
u{±L, z, t) — 0, the material is pulled in the left region 
{x < 0) and pushed in the right region {x > 0). When the 

u(H) 




FIG. 2: u{x, z, Ti/2) as a function of {x, z) in the case qm ~ 
0.08. 



deformation becomes sufficiently large, the shear stress 
CxziXjZjt) exceeds the yield stress (Jyo- Then, a plastic 
deformation occurs near the bottom (z = 0), as shown 
in Fig. O Due to this plastic deformation, the material 
remains pulled in the left region {x < 0) and pushed in 
the right region {x > 0) even after an external force a{t) 
becomes 0. Therefore, the normal stress a^xix, z^T^) is 
positive in the left region (x < 0) and negative in the 




FIG. 3: s{x, z, T3) as a function of {x, z) in the case aM — 
0.08. 



right region (x > 0), as shown in Fig. 0] In this manner, 



Oxx (1) 




FIG. 4; axx{x, z, T3) as a function of (x, z) in the case om = 
0.08. 

a material memorizes the effect of the force as a plas- 
tic deformation when an external force is applied to the 
material before drying. 

Next, we investigate the influence of a plastic defor- 
mation on the formation of drying cracks by comparing 
the time evolutions of {(Jxxix,t)) after starting the dry- 
ing process {t > T3) in the cases aM — and aM — 0.08, 
which arc illustrated in Fig. [51 and Fig. El In the case 
aM = 0, {(Txxix,t)) = just before drying (t = T3); 
however, in the case aM = 0.08, {(7xx{x,t)) is positive 
in the left region (a; < 0) and negative in the right re- 
gion [x > 0) by a plastic deformation. As the material is 
dried, the average normal stress {(Jxx{x,t)) increases in 
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a similar manner for both cases. However, until a crack 
is formed, the maximum values of (axxix,t)) in the case 
um — 0.08 are larger than those in the case aM = 0. 
Due to this large stress, a crack is formed earlier in the 
case aM — 0.08, than in the case um — 0- 
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FIG. 5; The average normal stress {axx 
om = as a function of x. 
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6: The average normal stress {axx{x,t)) in the case 
- 0.08 as a function of x. 



Here, we consider the correspondence between the re- 
sults of the viscoplastic model and the experiments. In 
the viscoplastic model, we consider a crack that is per- 
pendicular to an external force applied before drying. 
Due to the influence of a plastic deformation, the crack 
is formed earlier when an external force is applied to the 
material than when no force is applied. Therefore, we 
conjecture that in the experiments too, the perpendicu- 
lar crack emerges due to the influence of a plastic defor- 
mation when an external force is applied to the material. 



IV. QUANTITATIVE PREDICTIONS 

In the viscoplastic model, the position of the first per- 
pendicular crack deviates from the center of the material 
to the opposite direction of an external force applied be- 
fore drying as shown in Fig. |3 We wish to analytically 



express the position of the first perpendicular crack Xc 
as a function of the maximum value a^j of the external 
force. However, since it is difficult to study the partial 
differential equations ^ and (|2Jl, we simplify the equa- 
tions on some assumptions. 

First, we assume that an external force a{t) and a refer- 
ence strain rate c{t) vary more slowly than the relaxation 
of a displacement u{x,z,t) and a plastic strain s{x,z,t). 
Then, it is sufficient to calculate the static solution of 
the equations, given an external force Oi{t), a yield stress 
(7y(i), and a reference strain rate c(t). Second, we dis- 
cretizc the partial derivative of z as 



du{x,H,t) u{x, H,t) — u{x,Q,t) u{x,H,t) 



dz 



H 



H 



(12) 



da^(x,H,t) cj^^{x,H,t) ~ (Jx^{x,0,t) -<7x^{x,0,t) 



dz 



H 



H 



(13) 

where u{x,Q,t) = and (Jxz{x,H,t) = (see Eqs. ^ 
and (PJ). Based on these assumptions, the displacement 
at the surface U{x,t) = u{x,H,t), the normal stress at 
the surface Txx{x,t) = axx{x, H,t), and the shear stress 
at the bottom Txz{x,t) = (7xz{x,0,t) are determined by 
the equations 



dTxx{x,t) Txz{x,t) 



dx 



H 



+ a{t) = 0, 



Txx{x,t) = {X + 2^i) (^Ei^+c{t) 
\ ox 

Txz{x,t)=^,(^^^^-Six,t) 



(14) 

(15) 
(16) 



where S{x, t) is a plastic strain at the surface (s(x, H, t)). 
The time evolution of S{x, t) is as follows: if 



\Txz(x,t)\ < cry(i), 

S{x,t) does not change in time. However, if 
\Txzix,t)\ > avit), 

S{x, t) is determined by the condition 
\Txzix,t)\^aYit), 

which yields 

Uix,t) , arit) 
Six,t)^-^±—-. 



(17) 
(18) 
(19) 

(20) 



Here, the sign depends on the sign of Txz(x, t). By sub- 
stituting Eqs. ((TB|) and ((TB|l into Eq. lfTI|l . we obtain the 

equation of U{x, t) as 

(>^ + ^^^)^-^U+^S + a{t)^0. (21) 
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The functional forms of a{t), (Ty(t), and c{t) are given 
as Eqs. ®, ®, and iQ, respectively. From Eqs. © 
and the bomidary conditions are rewritten as 

U{L,t) ^U{-L,t) ^0 (22) 

for < i < T2 and 

r,,(x,t) = r,,(~L,i) = (23) 

for T2 < t. 

From Eq. the condition of a crack formation is 

rewritten as 



(24) 



Moreover, because T^xix^tc) has a maximum value at 
X — Xc, the equation 



d^xxiXc, ^c) 

9a; 



= 



(25) 



should be satisfied. 

From these equations, Xc can be calculated. The result 
is summarized below (see Appendix ^ for details of the 
calculation). First, we denote the threshold value of the 
external force by ayo, which is derived as 



ayo 



where 



ayo coshqL 
H {cosh qL - 1) 



(A + 2^)7?2- 



(26) 



(27) 



Then, if a^f < ayo, the position of the first perpendicu- 
lar crack Xc is expressed as 



Xc = 0. 



(28) 



In contrast, when ayo < um < ctYi, Xc is determined 



by 



coshgL — B{0,Xc)'^ 
D{L,Xs 



<JbqH 



ayo - aAiH 



sinh qXc 



sinh qL 
where ayi is defined as 

ayo{2 coshqL — B{L, Xg)) 



(1 — cosh qa;c cosh gi), (29) 



ayi 



(30) 



coshqL — B{L, Xg) 
Here, B{xi,X2) and D{xi,X2) are defined as 

B{xi,X2) — cosh(7(xi — X2) + qx2 sinh(xi — X2), (31) 

D{xi, X2) = sinh(7(xi — X2) + qx2 coshg'(a;i — 2:2). (32) 

Xs represents the region where a plastic deformation oc- 
curs as —Xs < X < Xg, which is determined by 

aMH+{ayo^aMH)B{L,Xs) = 0. (33) 



We further extract a simple expression of Xc by focus- 
ing on the region where a a/ is adjacent to ayo. Assuming 
that 



qXc < 1, 



qxs <. 1, 



(34) 



(35) 



and expanding Eqs. (|29|l and 1)33(1 to the first order of Xc 
and to the second order of Xg , we obtain 



(coshgL — 1)^ 
g^CTb cosh qL 



{aM - ayo)- 



(36) 



In the calculation of Xc, we simplified the original par- 
tial differential equations based on some assumptions. In 
order to confirm the qualitative accuracy of Eqs. (|28|l and 
((36|l . we numerically calculate the solutions of Xc for the 
original partial differential equations. In Fig. [7| and |H1 
we show the numerical solutions of Xc for the case where 
the parameter values are same as in Sec. IIIII (parameter 
set A) and for the case where the parameter values are 
set as H ^ 1.0, L = 15.0, A = 1.0, fi = 0.2, ayo = 0.03, 
CTb = 0.02, 7 = 1.0, B = 1.0, Ti = 300, T2 = 600, 
T3 = 500, and b = 0.000005 (parameter set B). There 
exists a threshold value ayo- Xc remains zero when the 
external force a^j is smaller than the threshold value ayo 
for both of the parameter sets. When the external force 
aM is larger than the threshold value ayo, Xc deviates 
from zero. In Fig. |51 we show the numerical solution of 
function of difference between the maximum value 
of an external force aM and the threshold value ayo for 
parameter set A. This figure clearly indicates Xc deviates 
from zero as a linear function of the difference between 
a A/ and ayo when aM is adjacent to ayo- These behav- 
ior qualitatively agree with Eqs. and though 
the values of ayo are different. 




0.045 0.05 



0.055 0.06 
ttM (X/H) 



0.065 0.07 



FIG. 7: The position of the first perpendicular crack Xc as 
a function of the maximum value of an external force um for 
parameter set A. 



Based on these results, we expect that in experiments 
too, the relation between the position of the first crack 
and the maximum value of an external force is expressed 
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FIG. 8: The position of the first perpendicular crack Xc as 
a function of the maximum value of an external force um for 
parameter set B. 
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FIG. 9: The position of the first perpendicular crack Xc as 
a function of difference between the maximum value of an 
external force um and the threshold value ayo for parameter 
set A. 



as Eqs. (|28|l and H3t)|) . The experimental confirmation 
of Eqs. (|28|l and H3t)|l supports our conjecture that the 
memory effect on the formation of drying cracks arises 
from a plastic deformation of the material. 

V. SUMMARY AND DISCUSSION 

In this paper, we model the formation of drying cracks 
in a viscoplastic material. In numerical experiments, we 
observe that when an external force is applied before dry- 
ing, a crack whose direction is perpendicular to the force 
emerges earlier than when no force is applied. This phe- 
nomenon occurs because of a plastic deformation. Based 
on this observation, we conjecture that a plastic defor- 
mation is the cause of the memory effect on the forma- 
tion of drying cracks. In order to check this theory, we 
quantitatively predict the position of the first perpen- 
dicular crack. If Eqs. (|28|) and (|36|l are experimentally 
confirmed, it may be concluded that a plastic deforma- 
tion causes the memory effect on the formation of drying 



cracks. 

Here, we remark on the neglect of some quantities 
in our model. A displacement w in the z-direction is 
not taken into account in our model. Moreover, normal 
stresses a^x, <^yy and azz sue not taken into account in 
the equation which represents the occurrence of a plastic 
deformation, i.e., Eq. ||2Jl. The reason for neglecting w 
is that w plays little role in the formation of the crack 
because a crack is caused by a normal stress axx and 
w hardly contributes to axx- The reason for neglecting 
f^xx , cFyy and azz is that we consider the situation where a 
plastic deformation occurs by a shear. To make sure that 
the neglect of these quantities does not affect our results, 
we simulated the model in which these quantities were 
taken into account. As far as we checked, we obtained 
qualitatively similar results with the model in present 
paper, for example, with regard to space distributions of 
a horizontal displacement u{x,z,t), stresses axxix, z,t) 
and axz{x, z, t) after applying an external force. 

In addition, we remark on the the results for another 
choice of the parameters and the expressions of otit), 
(Ty (t) and c{t). As far as we checked, we observe the same 
phenomenon that the first perpendicular crack occurs 
earlier and the position of the crack deviates from center 
for any choice of the parameter as shown in Fig. and 
Fig. |H| For any choice of the expressions of a(t), (Ty(i) 
and c(i), the same phenomenon are observed. Hence, we 
expect that our results are robust for the variation of the 
parameters and the expressions. 

Recently, Nakahara and Matsuo measured the rheolog- 
ical property of a powder- water mixture 0- From this 
measurement, they demonstrated that the mixture mem- 
orizes an external force before drying only when it be- 
haves as a viscoplastic material with a finite yield stress. 
Moreover, they report that only when the strength of 
the force is larger than a threshold, the mixture mem- 
orizes the force. The relation between the memory of 
the force and a finite yield stress supports our conjecture 
that a plastic deformation plays an important role in the 
memory effect on the formation of drying cracks. Fur- 
thermore, the existence of the threshold coincides with 
the result of our model. 
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APPENDIX A: THE ANALYTICAL SOLUTION 
OF THE POSITION OF THE FIRST CRACK 

In this appendix, we show the path to obtain Eqs. H28|) 
and (|29|l . which determine the position of the first per- 
pendicular crack Xc- This appendix consists of two sub- 
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sections. In the first subsection, we calculate S{x,Ti), 
and in the second subsection, we calculate Txx{x,t) (t > 
Ta). Substituting the expression of Txx{x,t) into Eqs. 
(|^ and we obtain Eqs. and In these 

calculations, we assume that Ti < T2 < T3 in order to 
simplify the calculations. 



1. The calculation of S{x,Ti) 

First, we evaluate the minimum external force ayo, 
which yields a plastic deformation. Assuming that a plas- 
tic deformation does not occur {S{x, t) = 0) and solving 
Eq. (|21() under the boundary condition (|22() . we obtain 



where 



U{x,t) = "^^^-^^ A _ coshg^ 
' ^ \ coshgL 



(A + 2^)i?2 



(Al) 



(A2) 



By substituting Eq. lAljl into Eq. we obtain a 

shear stress T^zix^t) as 



TUx,t)^a{t)H{l 



coshqx 



(A3) 



cosh qL ^ 

The assumption that S{x,t) = is valid if the equation 
T,z{x,t)<aYo (A4) 
is satisfied. This condition is equivalent to 

a{t) < ayo, (A5) 
I 



where 

ayo cosh qL 
" i7(coshgL- 1)' 

Hence, we find that if a^/ < ayo, 

S{x,Ti)=0. 



(A6) 



(A7) 



If ajvf > ctyo, a plastic deformation occurs. Then, 
assuming that a plastic deformation occurs in the region 
—Xs < X < Xs, we calculate S{x,Ti/2). From Eq. (jSUJ, 
t) is given by 



X < —Xs 

I U ^YO 

b[X,t) ^ < — -Xs < X < Xs 

x^ < X. 



(A8) 



By substituting this into Eq. H21() and noting Q!(ri/2) 
we obtain the equations of U{x, Ti/2) as 



(A 



^ 77^^ + = x< 



(A + 2^l) 



dx^ 



i/2 

gyp 

-p^C/ + aM — Xs < X. 



2M)^-4r + "M = -a;, < a: < X, (A9) 



^2 

By solving these equations under the boundary con- 
ditions H22(l and under the matching conditions that 
U{x,Ti/2), Txx{x,Ti/2), and S{x,Ti/2) are continuous 
at a; = ixs, we obtain the expression of U{x,Ti/2) 
and the equation to determine Xs- The expression of 
U{x,Ti/2) is 



U{x,T^/2) = { 



1 B(x, -Xs) 

AsH o, T o, o'yoH 



2/2 2\ I 

2^^ -•"^)+ M 

aMH^ 



AsH 



Xg 

Xs < X, 



(AlO) 



where 



and 



As = cryo - ttM-ff 



(All) 



B(x\,X2) — coshg(xi — X2) -f (7a;2 sinh(xi — x-i). (A12) 
The equation to determine Xs is 

aMii + AsB(L,Xs)=Q. (A13) 



Substituting Eq. IjAlOp into Eq. IjASp . we obtain 
S'(x,Ti/2) as 







A 



S{x,T^l2)={ ^q\ 



J2./2 •2\ 



X X r. } X i 



< X < Xs (A14) 







Finally, S{x,Ti) is calculated. First, we calculate 
U{x,Ti) on the assumption that 



Six,Ti) = S{x,T,/2). 



(A15) 
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Substituting Eq. l|IT5|) into Eq. (ETJ and noting a{t) = 0, we obtain the equation of U{x,Ti) as 
I 



(A + 2/i)|^ - -^U + ^q^{x^ -xD^Q ~Xs<x< Xs 



A, 



(A + 2/i) 



dx"^ 



2H 



(A16) 



u = 



r 



By solving this equation under the boundary conditions and the matching conditions, we obtain U{x,Ti) as 
I 



U{x,Ti) = { 



^(Bix.-x.)-Bmx.f-^ 
/i \ cosh qL 

AsH / ^ q^jx"^ -xl) _ ^ ^ . coshgx 

/i \ 2 ' ^ coshqi 

/ , . , coshqx' 
B{x, Xs) - B{L, Xs) 



cosh qL 



Xs < X. 



(A17) 



r 



Substituting this equation and Eq. (|A15p into Eq. (fl^ . we get Txz{x,Ti) as 



-B(x, -Xs) - B{L, Xs) 



Txz{.x,Ti) = < 



cosh qx 

cosh qL 
. ( „, N coshoa; 

. / „ , s „ , N cosh qx 

As B(x,Xs) - B(L,Xs) — 

V ^coshgL 



Xs X 3 



r 



(A18) 



The assumption that S{x, Ti) — S{x, Ti/2) is vahd when 
the equation 



|21rz(x,Ti)| < avQ 



(A19) 



is satisfied, because S{x,Ti) — S{x,Ti/2) imphes that a 
plastic deformation does not occur for Ti/2 < t < Ti. 
This condition is equivalent to 



where 



an 



cryo(2 coshqL — B{L, Xs)) 



coshqL — B{L, Xs) 
Hence, if 

avo < CUM < CHYl, 

then S{x, Ti) is expressed by Eq. HA15|I . 



(A20) 
(A21) 
(A22) 



2. The calculation of T^x(x,t) for t > T3 

In this subsection, we calculate T^xix^t) for i > T3 in 
the cases um < ceyo and ayo < ctM < otyi- 

If ttM < oiYo, then S{x,Ti) — 0. Since cry — 00, 
5(2;, t) = for t > T3. Substituting S{x,t) = and 
a{t) — into Eq. (|^ . we obtain the equation of J7(a;, i) 
for i > T3 as 



(A + 2^) 



f7 = 0. 



(A23) 



By solving this equation under the boundary conditions 
we obtain U{x, t) for i > Ts as 



U{x,t) 



c{t) sinh 
(7 cosh qL 



(A24) 
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case aM < Oiyo as 



Substituting this into Eq. (|15|l . we get Txx{x,t) in the for t > T3, where S{x,Ti) is expressed by Eq. ljA15|l . 

Substituting Eq. (|A26II into Eq. (EU and noting a{t) = 
0, we obtain the same equation of U{x,t) as Eq. ljA16p . 
By solving this equation under the boundary conditions 
(|23|l and the matching conditions, we obtain U{x,t) for 
t > T3 as 



T,x{x,t) = {\ + 2^l)c{t) 1 



coshqx 
coshqi 



(A25) 



If ayo < aM < otYi, S{x,Ti) is expressed by Eq. 
(|A15|) . Since ay = 00, 



S{x, t) — S(x, Ti) 



(A26) 



U{x,t) = < 



c{t) sinh qx 




q cosh qL 




c{t) sinh qx 


AsH 


q cosh qL 




c{t) sinliga; 


AsH 


q cosh qL 





B{x, -Xs) 
q^ix^ 



coshqx 
D{L,Xs) . 

smh qL 

coshqx 

D[L,Xs) . 

smh qL 



cosh qx 
B[x,Xs) - D[L,Xs)^- — - 
smh qL 



Xg X Xg 



(A27) 



r 



where 

D{xi,X2) 



smhq{xi — 2^2) + 92^2 coshg(a;i —2:2). (A28) 



Substituting this into Eq. H15() . we get Txx{x, t) for i > T3 
in the case ayo < ctM < ayi as 



(A + 2^l)c{t) i 1 - 
(A + 2M)c(t) (1 - 
(A + 2^)c(t) (1 



cosh \ 



cosh qL 
cosh gx 

cosh qL 
cosh gx 



+ 



qH 
A, 



coshgL / qH 



sinh qx 
D^x, -Xs) - D{L, Xs) — 



qx 
D{x, X 



sinhqx 
D[L,Xs)- 



sinh qL 



sinh qL 

sinhqx 
D{L,Xs) . 

smh qL 



X < —Xs 
~Xg X s 



(A29) 
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